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Nomenclature 


A..( r .t): 

1J ** 

second-order anisotropic tensor 


two-point time-dependent second-order velocity correlation 


tensor 

c „ : 

decay constant defined in Equ. (12) 

F,C**),F(—): 

polynomial basis functions for Q(***) 

1. ( r ,t): 
ij~ ’ 

second-order isotropic tensor 

M ,M : 
1 2 

mean-flow wind shear gradients 

p , n (q,r): 

homogeneous polynomial of order m in the variables, q,r 

Q(q,r,t): 

third-order three-point scalar correlation function 

R(r,t): 

normalized scalar velocity correlation function 

S.. (r.t): 

ijm v ~ 

third-order two-point velocity correlation tensor 

U : 

i 

total flow velocity 

U: 

i 

mean-flow velocity 

V: 

airplane flight speed 

a,b: 

time-dependent parameters characterizing anisotropic 


turbulence 

f(r,t)/g(r,t): 

longitudinal/transverse correlation functions for i 


turbulence 

k(r,t): 

scalar function peculiar to third-order velocity moments of 


turbulence 

X : 

separation vector 

r: 

magnitude of separation vector 

u.: 

1 

turbulent flow velocity 

w(x,t): 

turbulence downwash 

x,z: 

Cartesian coordinate variables 

A(ty. 

integral scale for isotropic turbulence 
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AA(i): 

change in integral scale due to anisotropy 

A{k, t): 

energy-scaling factor 

a: 

airplane glide-slope angle 

P,Y- 

dimensionless constants 

k: 

non-dimensionalized wave-number 

i: 

imaginary element 

<r(t): 

isotropic turbulence intensity 

A<j(t): 

change in intensity due to anisotropy 

oj: 

angular frequency 

0: 

dummy variable 

<t>{K, t): 

power-spectral density 

<f> a \ic, t): 

power-spectral density for isotropic turbulence 

4> <a W): 

anisotropic part of power-spectral density 

r: 

time-lag 

n- 

generic velocity scale factor 

*>• 

inverse length scale factor 
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1. Background 

Meteorological events such as thunderstorms and unstable frontal systems have 
long been considered dangerous from the perspective of aviation safety. 

Investigations into the alarming number of recent low-level aircraft accidents 
involving thunderstorms around airports have revealed one meteorological 
culprit — a strong blast of air directed toward the ground which has come to 

be known as a "microburst". 

The mechanism by which the microburst is created is a complex and sometimes 
violent one. From the investigation of a particular airliner crash at JFK 
airport, much was learned about the convective atmospheric dynamics which 
cause the microburst. It was postulated by Fujita* that a violent atmospheric 

disturbance may occur when moist upper air drops precipitation through a 
relatively dry air layer below it. As the precipitation evaporates, the dry 

layer cools: consequently, a stream of dense cold air suddenly replaces the 

ground-heated low-altitude air. Furthermore, as the resulting downdraft 
penetrates through the existing horizontal flow, turbulence is invariably 

generated. The downdraft eventually impacts the ground and spreads out 

radially. During the evolution of a microburst, the dimensions of the 

low-altitude phase may be on the order of 15-30 meters high with a radial 

spread of 2.5 kilometers. Realization of its relatively small meteorological 

size serves as the motivation for the term "microburst". Wind velocity 
information subsequently obtained from flight recorders has sustained Fujita's 

explanation. The typically distinguishing characteristics of the microburst 
are winds which are localized and variable. Abrupt changes in wind speed 
and/or direction over a short atmospheric distance are referred to in the 
literature as "wind shear". For purposes of this study, however, wind shear 

will be defined as a flow which has a spatially nonconstant mean-velocity 
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profile. Flight through a "wind shear-infected area" will tax both pilot and 

aircraft since related velocity gradients may drop airspeed to critically low 
levels. In summary, microbursts and wind shear pose an immediate threat to 
flight operations at low altitudes due to their adverse effects on aerodynamic 
lift/aircraft response. 

The scenario of microburst encounter by an aircraft in the landing mode is 

shown in Figure 1. First the aircraft senses an uncharacteristically strong 
head wind. Although turbulence typically accompanies the increased airspeed, 
pilot confidence is high during this phase of the encounter since aircraft 
performance increases with the additional aerodynamic lift. Because the 

intent is to land the aircraft, the pilot intuitively trims. Eventually, the 

head wind becomes a tail wind, and at this point performance is immediately 
and seriously degraded since the distinct head-to-tail wind swing has reduced 

airspeed enough to possibly lose flying speed. In the last ten years, an 

estimated 575 people have been killed in commercial aircraft accidents in the 

United States alone due to the microburst phenomenon. Likewise, 
microburst-induced accidents have also been reported in Great Britain, 
Germany, France, Italy, Australia, and Japan. Consequently, a fundamental 

understanding of microbursts, related wind shear and turbulence, is a 

high-priority research issue. 

Flight science authorities agree that the "solution" to the wind shear 

"problem" is multi-faceted. The existing research programs of NASA and FAA 
focus primarily on three elements: 1) hazard characterization: 2) sensor 

technology; and 3) flight management and operations. Hazard characterization 
is the study of the physics of the microburst phenomenon. Inherent within 
this phase of the research effort is analysis of aircraft aerodynamics in wind 
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shears and in heavy rain. Wind shear velocity profiles, rainfall effects, and 
turbulence models are important contributions from this component of the 
programs. Sensor technology deals with the prediction and detection of 
potentially dangerous meteorological events; it further subdivides into 

airport-fixed and airborne sensor technologies. Next Generation Radar fields 
(NEXRAD) are presently being investigated by NOAA, FAA and USAF for future use 
as airport radars. The other detection device is the in situ look-ahead 
sensor. This onboard computer system is being devised to scan the forward 

atmospheric environment, evaluate a "hazard index", and annunciate a warning 
to the pilot in a time period adequate to ensure either avoidance or escape of 
the threat. Look-ahead sensor technology is beyond state-of-the art, and is 
not expected to be operational until the mid-1990's. The third component of 
microburst research is flight management. Simply put, this phase concentrates 
on informing the pilot on how to safely get the aircraft out of a microburst 
encounter subject to aircraft performance constraints. An important element 
of such studies is flight simulation. Simulators allow pilots to experience 

threatening wind shears in a controlled environment, with intent to better 

prepare them for potential real-life events. Realistic wind shear 

representations in flight simulators are practical and economical, and can 
help crews coordinate their "escape" efforts in critical situations. Realism 

in the simulation of flight through hazardous atmospheric environments has 
improved with the introduction of lateral and vertical winds, as well as 
vortex and turbulence influences. The focus of the present work is in the 

hazard characterization aspect of the programs; its principle objectives are 
to model the fluid-dynamical characteristics of turbulence associated with 

low-level sheared mean-flow, and quantify the effect of these characteristics 

on aircraft response. All of these issues will be herein addressed from the 

perspective of p ower s pectral density (psd). 


6 


2. Introduction 

Since microbursts pose a hazardous threat to flight operations only at low 
altitudes, the role that the ground plays as a fixed boundary is therefore 

important in terms of turbulence modelling. The presence of the ground 

generates in the microburst a brand of turbulence which when sensed by the 
aircraft during its take-off or landing is both nonstation arv and a nisotro pic. 
The objectives of this investigation are thus to model the statistical 
characteristics of same, and gain a quantitative understanding, in terms of 
scaling effects, of the stochastic nature of the turbulence sensed by an 

aircraft during a microburst encounter. 

The task of mathematically describing the turbulence associated with a typical 
microburst (shown in Figure 1) is straightforward. The pioneering works of 

A A J 

Taylor , Karman and Howarth , and Robertson^ have illustrated the well-defined 
fashion in which the concepts of two-point velocity correlations can be used 
to quantify the kinematical aspects of isotropic turbulence; Batchelor^ 

extended these concepts to the case of axisvmmetric turbulence that which 

has a defined form of anisotropy. Axisymmetric turbulence, with its tendency 
to prefer one spatial direction, is the logical candidate for modelling wind 

shear turbulence (with its particular anisotropy resulting from the microburst 

impacting the ground). Although classical axisymmetric theory produces 

turbulence velocity correlation functions that are position-independent 

(homogeneous) and time-varying, the turbulence under consideration is assumed 
inhomogeneous and in a steady-state ("frozen") condition; nonetheless, it can 
be easily shown** that correlation functions for frozen, inhomogeneous 

turbulence may be converted into time-evolution equations for homogeneous 

turbulence by duly scaling out a characteristic velocity. Accordingly, 
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C..(x . r ) 

r <v <v» 




physics of the inhomogeneous, steady-state turbulence aligns with the physics 
of classical turbulence; the ( r ,t)-notation rather than the (t, r )-notation is 
the accepted standard. 


For the specific case of an airplane in either the take-off or landing mode, 

n 

recent work , which incorporates axisymmetric theory, has resulted in the 
"airplane-sensed" two-point velocity autocorrelation for anisotropic 
turbulence downwash, viz. 


C 33 (r,t) = <w(x,t) w(x + r ,t)> * {a g} + {a|r| + b} ; 


( 1 ) 


here the magnitude of a ~ 0(a), where a is the airplane glide-slope angle, b 
is a parameter (to be defined later) which can be either positive or negative, 
and g is the transverse correlation function of isotropic turbulence. 

Equation (1) is based on the assumption that during take-off or landing the 

aircraft senses only a small amount of anisotropy in the (necessarily) large 
amount present in the microburst. Note that the velocity correlation 
separates into isotropic and anisotropic components, and takes the general 

form, C„<r,,) ~ {jSL} + 


/isotropic 1 _ 2 

(correlation/ g 


and 


- ’M + »• 


The contributions of isotropy and anisotropy to the total correlation are 
shown for b > 0 only in Figure 2; recall that b can be either positive or 
negative. 


An important aspect of the ensuing analysis is how turbulence "scales" in a 
inicroburst. A common correlation used for nonstationary processes is the 
uniformly modulated, viz. C(r,t) = a 2 (t)R(r) where cr 2 (t) is the time-varying 

modulator; implicit in this correlation is the constraint imposed on the 
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g 

integral scale — that of time-invariance. In recent studies , however, a more 

general form has been formulated which incorporates time-dependence in both 
intensity and integral scale, viz. C(r,t) = tjR(£r) where rj = //( t) as a 2 ( t) and 
4 — £(t) as 17/1 (t). Allowing the correlation to evolve in such a manner 

permits its geometric shape to distort in time, while its r-algebraic form is 

retained. Figure 3 illustrates the "self-preserving" behavior of such a 

correlation function. With C(0,t) in this case * a (t), intensity for 

self-similarity is clearly an ordinate-axis scaling effect while the 

normalized correlation function, R(£r), automatically incorporates the 
integral scale/abscissa-axis scaling effect. Although self-preservation is 
introduced here strictly from the standpoint of "scaling-effects", its 

hypothetical existence and subsequent analysis can be be justified on purely 

"dynamical system" grounds (see Appendix). 

3. Power Spectral Density Analysis 

The power spectral density is a common way of illustrating energy distribution 
versus frequency, and is an immensely popular tool in both gust analyses and 
turbulence simulations; mathematically, it is the Fourier transform of the 

autocorrelation function. The psd for the aforementioned down wash formulates 
as 

0(<y/V,t) = J <w(x.t)w(x + r,t)>e' /Wr/v dr = J {a 2 g + a|r| -I- b}e 4a>r/v dr, (2) 

where i 2 = -1. Since the correlation function consists of isotropic and 

anisotropic components, so too will its transform; specifically 0(K,t) = 

0 a> (K,t) + 0 ,a> (K,t), where 

0 ( V.t) ■ <t 2 A and 0 <8) (/c, t) = 4a 2 A (^41 

(1 +rc ) 1 J (Kb/a/1) 2 

Here <f>(K, t) is the total psd sensed by the aircraft as it flies through the 
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microburst, k = coA/V a non-dimensionalized wave number, a the isotropic 
intensity, A the longitudinal correlation length (integral scale), and AA the 
change in integral scale due to anisotropy. The scaling effects discussed 
earlier are evident upon inspection of Equ. (3). 

The distinct advantage that the psd representation has over its transform 
partner, the autocorrelation function, is its more defined statement of the 
modelling parameters of anisotropy, viz. AA, a, and b, in relation to total 

turbulent energy. For a glide-slope which is maintained during a microburst 
encounter, it follows that a * constant since (recall) magnitude {a(t)} ~ 

Q(«). Hence, b plays the more important role in terms of modelling the 

anisotropy due to the boundary presence. Since b is related to intensity by 
a + b = (ct + A a) , Ao ~ increment in intensity contributed by anisotropy, for 
neglected higher-order terms b becomes 

» = * 2 {H 2 - + [#“)} - 2 ° 2 [¥) = m*’ ( 4 ) 

note that for pure isotropic turbulence, b ® 0. One way to quantify the 

collective effects of AA, a, and b on the airplane-sensed turbulence is to 

introduce an "energy scaling factor", defined as A{K,t) = 

To analyze the effect of AA/A on A(K,i), set b/(a/l) <* constant; typical 
values of a = 5 ft/s, A = 1000 ft, V = 225 ft/s and a ~ Q(a) = Q(3°) = Q(0.052 
rad) are chosen, and Aala is allowed to vary up to 0.5 (a feasible value for 
wind-shear turbulence severity). Incorporating Equ. (4) gives b/(a/l) - 

Q(0.5). Figure 4 shows that for AA/A = 0.1, the energy scaling factor is as 
high as 35% for 0 < k < 2 (0 s co < 0.45 rad/sec) — a frequency range which 
falls within that specified for the guidance and control phases of an aircraft 
mission involving a wind shear. Hence, a realistic estimate of the energy 
transferred to the airplane requires that anisotropy be accounted for in both 
turbulence models and related flight simulations. 
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4. Constraints Imposed by Navicr-Stokes 

Knowledge of the sheared mean-wind/turbulence interplay requires that the 
governing equation of the fluid be considered. The Navier-Stokes equation for 
a turbulent, high Reynolds number, steady-state fluid with both body and 
pressure forces neglected is 


d(U U ) 

1 m 

cFx 


~ 0. 


With U. == U + u , the two-point correlation method^ yields 


. au . au ' 

1- C . + a— - C - 2 Jr- S + — i- c =0 

01 it OX mi Or imj / lm 

in m OX 


(5) 


( 6 ) 


where C C (r.t) = <uu'>, S <^S(r.t)=<uuu'> and (cf . Figure 

ij ij ~ i j imj linj ~ l in j 

1) r = (r, 0, a r); index contraction results in 

au. au 


4c. + C . + t C. - 2 -J— S. . = 0. 

Ot li OX nil o _ lm Or 11m 

in OX in 


(7) 


Equation (7) requires the solution of separate isotropic and anisotropic 
dynamical problems, viz. 

a a 

( 8 ) 


"lr 1 ii ' 2 ~lr S imi ~ 0 


and 


au. 


au. 


T A u + 


ax 


c . + 


ax 


-4c , = 0 , 


(9) 


7.10 


I + A = C . The basis for a scaling comparison (the isotropic 

ij ij ij 


where 

case) will be considered first 


Introducing the f and g functional representations for the correlation tensors 

gives a result whose first integral is the Karman-Howarth equation for high 

o 

Reynolds number isotropic turbulence , viz. 

a, «*> - (>— + 4 r) “ °- < 10 > 

where k = k(r,t) is the scalar function required to specify third-order 
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velocity product moments. Since turbulence self-interaction is not altered by 

self-preservation, it follows that k(r,t) also has such an evolution with 
time, i.e. k(r,t) k(£r) = k(r//l). Incorporation of this feature into the 
r-integration of the Karman-Howarth equation results in 

-jf = <7 5 C k . (II) 

where 

C fc * J 4. d0 = negative scale-independent constant, (12) 

and 0 = r/A; a specific numerical value of C fc will be determined later. 
Solution of Equ. (11) gives the time-dependence, 

A{t) = J cr 3 (t)dt, (13) 

a 2 {t) 

for the integral scale. The power-law form for isotropic intensity, cr 2 (t) ** 

t n , n = ±1, +2, ... , converts Equ. (13) into 

2 C 

Ait) - t (n+2)/2 , N = ±1, ±2, ... ; (14) 

for n = - 1 , Ait) ~ 0(t 1/2 ). 

To consider the anisotropic part of the problem, viz. Equ. (9), gradients must 
first be defined. The two wind-shear profiles chosen for analysis are 
illustrated in Figures 5 and 6, and will hereinafter be referred to as "Shear 
Flow I" and "Shear Flow II", respectively. Note that both profiles have a 
"head-to-tail wind swing", a phenomenon frequently observed in microburst 
encounters. The mean velocity for Shear-Flow I is U. = {U(z),0,0} = {MjZAO} 
as shown in Figure 5; Mj is negative on the headwind side of the microburst 
and positive on the tailwind side. Analysis of Shear Flow II is restricted to 
the two-dimensional case with U. = {U(x),0,W(z)} and gradients of the form 
- = M 2 > 0. Incorporation of Shear-Flow I characteristics into 
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Equ. (9) yields 


(a j r | +b) - MaV + M^ra 2 - 0 ; 

integration in r from 0 to oo produces 

4 r {o\t)AA(i)} - M^Oo/Kt) = 0 , 

oo 

where J (a|r| + b) dr = a 2 J/l(t). The solution to Equ. (16) is 
o 

M a 

AA( t) = — | cr 2 (t)/l(t)dt. 

(7 (t) 

For cr 2 (t) t” and the integral scale of Equ. (14), Equ. (17) becomes 


AA(l) 


4M C, a 

1 k 


(3n __ + 2]T3n 4- 4) 


t (n+4)/2 , i I 9 

t , n = ±1, ±2, .... 


(15) 


(16) 


(17) 


(18) 


specifically, AA{t) ~ Q(t 3/2 ) for cr 2 (t) ~ t' 1 . A similar analysis for 

Shear-Flow II yields 


AA( t) 


4M 2 C t 


nn'+T)TTrm) 1 


(n+4)/2 „ _ 


, n = ±1, ±2, ... 


Like Shear Flow I, Shear Flow II has AA(t) ~ Q(t 3/2 ) for <r 2 (t) ~ t 1 . 


5. Determination of C, 

k 

The theory of turbulence strictly prohibits the presence in the analysis of 
any adjustable constants; therefore a full closure of the dynamical equation 
is not achieved until a value for C is determined. The fact that C is 

k -K 

formed from an integral of k(***) (cf. Equ. (12)) explicitly suggests that 
aside from some fundamental limiting constraints, viz. i-4 & (k(«»*)} ~ Q(rj* 
and & (k(*»»)} = 0, the exact functional form of k(»**) is not critical, 
i.e. k(**») itself is not necessarily a unique function. All that is required 
for any and all possible k(***), say k(***) k N (***), N = 1,2,3, . . ., is 

oo 

that (4 j r' l k (r,t)dr) = C be independent of N. A preliminary step is 
o N k 

therefore to develop a representative algebraic form for k(r,t) ~ k(r//l). 
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To accomplish this while simultaneously saying as little as possible about 
k(»»») itself, it is here proposed to kinematically describe the three-point 
correlation, <u(x,t)u(x + q,t)u(x + r,t)> = Q(q,r,t), where u(x,t) is only one 
component of the three-component turbulence velocity field and x is a 
one-dimensional position vector, i.e. x = (x,0,0), etc. Note that for 

self-preservation, this correlation accordingly scales as Q(q,r,t) = 
<7 3 (t)s(qM,r/y1), where s(0,rM) m k(r//l); furthermore, the functional form of 
Q(** # ) is such that it must satisfy certain symmetry properties, viz. Q(q,r,t) 
= Q(r,q,t) = Q(-q,r-q,t)**. Since Q(***) expanded in a Taylor's series about 
q = r = 0 produces 

Q(q,r,t) “ a 0 + a ( q + b^ + a 2 q 2 + c 2 qr + b 2 r 2 + ... , 
it follows from the required symmetry that the time-dependent parameters, a , 

in 

b , and c , m = 0,1,2,..., cannot be completely arbitrary. In fact they must 

1 2 

each be such that eventually Q(***) expresses as 

Q(q,r,t) = E A m (t)P (q,r) 

m - 0 

where each P (•••) is a homogeneous polynomial in q, r and individually must 
themselves satisfy the same symmetry as Q(q,r,t), i.e. P (q,r) = P (r.q) = 

m in 

P m (-q,r-q), m = 0,1,2,... . Fortunately, there are only two (2) fundamental 

such forms, viz. F ( (q,r) = q 2 - qr + r 2 and F 2 (q,r) = 2q 3 - 3q 2 r - 3qr 2 + 2r 3 , 

since any other higher-order such form can be represented as a combination of 
Fjfy**) and/or F (•••); there is no linear homogeneous algebraic form in q,r 
which satisfies the required symmetry. For example, the fourth-order 
polynomial is identical to (F (q,r)} 2 , while the fifth-order one is identical 

to (F i (q,r)F 2 (q,r)}. The existence of such a complete set of polynomial basis 
functions for any sequence of systematically generated algebraic forms was 

1 'J 

originally established in a theorem formulated by David Hilbert in 1890 . 

The theorem need not be quoted here, but the essence of it is that if an 
infinite sequence of algebraic forms (the sequence here is P Q (q,r) = 1, 
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J 

Pj(q*r) = o, P 2 (q,r) = q 2 - qr 4- r 2 , P 3 (q,r) = 2q 3 - 3q 2 r - 3qr 2 + 2r 3 , etc.) 

is constructed according to some general rule (the "rule" here is the required 

q,r-symmetry), there comes a point in the construction where the set of all 

available independent such forms is exhausted and any form constructed 

thereinafter contains absolutely no new information, being merely a 

modification/repetition of what has already preceded. What this means for the 

case at hand is that the most general sixth-order form must be a linear 

combination of the products, F f (q,r) and F 2 (q,r), while a seventh-order such 

form is {F j 2 (q,r)F 2 (q,r)} > A simple appropriate algebraic expression for 

14 

Q(q,r,t) which for q s 0 fits data reported by Stewart is thus 

, AF (q,r) 

Q(q,r,t) « <7 (t) 

{1 + BF^q.r)} 2 

where A = fiA 3 (t), B = yA~ 2 {t), and fi,y ** dimensionless constants; 

accordingly, 

k(r,t) ~ k(r//l) - WWA) 3 

{1 + y(r!Af} 2 

The specific fi and y for the data sets highlighted by Stewart (loc. cit.) are: 

Curve 1 (cf. Figure 7); = -1.48, y — 7; Re = 5300; 

Curve 2 (cf. Figure 8); /? = -3.0 , y = 12; Re = 21,200; 

Curve 3 (cf. Figure 9); fi — -7.42, y = 22; Re = 42,400; 

The value of C fc which most closely approximates all three cases is C fc - 

-0.167T. 

6. Conclusions and Concluding Remarks 

Since integral scale provides a measure of turbulence randomness^’ 
realistic estimates of this scale are essential to ensure related flight 

simulation fidelity. The role played by anisotropy in a wind shear encounter 

is quantified in the correlation length, L ■ A(t) + AA{t). As an aircraft in 

the landing mode penetrates the headwind portion of Shear Flow I, AA > 0 
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since M and « are negative in this phase of flight. The turbulence sensed 
there by the aircraft is therefore less random than that predicted by 
isotropic theory. The tailwind phase of such an encounter provides AA < 0 — a 
more random turbulence in comparison with isotropy. A similar landing mode 
analysis of Shear-Flow II yields negative A A on both sides of the microburst, 
thus making the attendant turbulence more random than the isotropic case in 
all phases of the encounter. 

The significance of anisotropy in microburst turbulence models has herein been 
considered from the perspective of self-similarity in the turbulence power 
spectral density. The considered energy scaling factor clearly defines the 
inadequacy of using isotropic theory only to estimate the necessarily 
anisotropic turbulent energy in a typical wind shear. The formulated psd 
representation of the turbulent downwash is a major step towards enhanced 
simulation of low-level turbulence, with the relative importance of the 
characteristics considered above becoming fully quantified only upon eventual 
flight-simulator implementation. 
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Appendix: The Case for Self-Preservation 


Nonstationary behavior whose time-evolution can be effected by two 

time-dependent scaling effects is important not only because of its relatively 
simple features, but also because of its immediate practical application. For 
nonstationary processes the two-time correlation function, Q(t ( .t 2 ) = 

<w(t j )w(t 2 )> , can be considered as a two degree-of- freedom system in which t 
= (t ( + t 2 )/2 and t = t - tj play the role of uncoupled principal 

coordinates; indeed in the (t,r)-coordinate reference frame the two 

degree-of-freedom system at hand is instead duly converted into one whose 
physical characteristics accordingly vary with time, i.e. Q(t ( ,t 2 ) Q(t,r) -» 
C(t,r/V) C(t,r), where V rr. And even though time-varying systems can be 
regarded as filters which "select" eigenfunctions in the pure sense rather 
than just simple sinusoids, the natural tendency in their analysis is to treat 

them as systems whose properties do not change "appreciably" with repect to 
time. The underlying motivation for this approach is that classical modal 
analysis is not explicitly restricted to systems which are time-invariant 
only; it is only necessary that the differential equation governing the 
system, as well as the boundary conditions, be separable into a function of 

time alone and in this case a function of non-dimensionalized length alone, 

i.e. that temporal changes in the system be in some sense "uniform". The 

suggested separation for dynamical systems is thus the essence of 
self-preservation for nonstationary behavior in random processes. The 
mathematics of same proceeds as follows. Consider the most general case of 
time-varying intensity and integral scale, viz. C(t,r)=cr(t-T/2)cr(t + T/2)R(t,r), 

r = r/V. If is expanded in Taylor’s series about r = 0, the result for 

C(t,r) is 

C(t,r) - cr 2 (t) { 1 - T 2 /4T 2 }R(t,r) 

where only linear terms in r have been retained for a( t ± t/2), and T = 
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! ~ 

a(t)[^] is the time-dependent scale over which a(t) may be assumed to vary 
linearly with r, viz. 

cr(t ± t/ 2) “ a(t){l ± r/2T}. 

Since R(***) scales with A(t), i.e. R(t,r) R(tM) 1 , and scales with 

T, the modified expression, 

C(t,r) - a 2 (t){l - (r/2/l) 2 (/l/T) 2 }R(** # ) , 

allows for analysis of the coupling between the two scales, T and A. Note 
that as far as R(***) is concerned, R(***) “ 0 for t > A, so that the major 
contribution to C(t,r) from R(»**) itself is from values of r < A. Therefore, 
if T » A, it follows that C(t,r) - <7 2 (t)R(***), i.e. the r-dependence of 
C(t,r) is determined by R(***) only and is unaffected by the r-dependence of 
{o ( t - r/2)cr(t 4- r/2)}. 


f A 


<v A (t) in this expression considered as a "time-length" scale. 
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the explicit two scaling-effects feat 



of anisotropy, i.e. AA/A. 
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Figure 7: Triple— velocity curve, k( . . . )» obtained by Stewart; broken line is 




experimental curve 
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Figure 8: Triple-velocity curve, k(. . .)> obtained by Stewart; broken line is 




experimental curve 
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Figure 9: Triple-velocity curve, k(. . .), obtained by Stewart; broken line is 
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